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Моделирование динамики и балансировка гибкого ротора по л-й форме 
изгиба оси 


О. О. Полушкин, О. А. Полушкин 
(Донской государственный технический университет) 


Посредством моделирования динамики гибкого ротора установлены закономерности, определяющие упругий 
прогиб его оси по собственным формам изгиба и связанные с этим распределённые дисбалансы. Разработа- 
на оригинальная методика идентификации параметров модели и балансировки таких роторов по собствен- 
ным формам изгиба. 
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Введение. Вследствие конечности изгибной жёсткости оси любого реального ротационного агре- 
гата (ротора) машины он имеет бесконечное множество критических скоростей и» вращения, 
соответствующих каждой из / = 1, 2, ..., П, ... собственных форм изгиба его оси [1]. Там же обос- 
новывается, что к гибким следует относить роторы, максимальная эксплуатационная скорость в, 
которых отвечает условию ‹, > 0,8571. Гибкий ротор балансируется по л-й форме изгиба [2], 
если р, ) <, <. Прил=1%,, )= 0 и гибкий ротор балансируется по 1-й форме изгиба, 


п- (п- 


если 0,8574, „ <, <@„„. 


кр1 


У каждого реального ротора, сбалансированного на низкой скорости в; < и, по тради- 


Кр1 
ционной технологии динамической балансировки с коррекцией дисбалансов в двух плоскостях 
коррекции у опор, объективно наличие исходного неупругого искривления оси, названного в [1] 
линией эксцентриситетов масс ротора. Эта линия представляет собой пространственную кривую 
отклонений (эксцентриситетов) центров масс элементарных сечений ротора, перпендикулярных 
его оси, от оси его вращения (прямой, соединяющей центры опор). Эту кривую можно разложить 
по собственным формам исходного неупругого изгиба [3], каждая из которых лежит в некоторой 
плоскости, содержащей ось вращения ротора. 

Вследствие исходного изгиба оси ротора по каждой из собственных форм дисбалансы от 
него распределяются неравномерно по этой форме. С приведением ротора во вращение со скоро- 
стью >, усилия от этих дисбалансов ведут к упругой деформации оси по всем собственным 
формам её изгиба. Чем ближе @ к @ни /= 1,2, ..., П, ..., тем значительнее эта деформация по Ёй 
собственной форме. При & =@&„„ наступает резонанс по Ай форме изгиба. 


Существующие методы балансировки гибких роторов (высокочастотная балансировка) не 
имеют общей теоретической основы. Учёт упругих изгибных деформаций оси гибкого ротора при 
его балансировке и в эксплуатации «крайне усложняет процесс балансировки и затрудняет уста- 
новление общей методики балансировки, пригодной для гибких роторов любой категории» [3]. 

Представленные ниже решения позволяют устранить эти недостатки и создать общий ал- 
горитмический подход к последовательной балансировке гибкого ротора по каждой из / = 1, 2, ..., 
п собственных форм изгиба. 

Моделирование динамики ротора, сбалансированного по / = 1, 2, ..., (л — 1)-й собст- 
венным формам. Объект исследования — двухопорный ротор в виде однородного гладкого кру- 
гового цилиндра с исходно изогнутой осью по всем собственным формам. Этот ротор предвари- 
тельно сбалансирован динамически на низкой скорости вращения 4, <, по традиционной 
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технологии в двух плоскостях коррекции у опор. Кроме того, он последовательно сбалансирован 
по /= 1,2, ..., (л — 1)-й собственным формам по специальным методикам, о которых речь пойдёт 
ниже. Полагаем проведённые балансировки идеальными, поэтому вращение такого ротора со 


скоростью &„„,) <@<а,„ не вызовет дополнительных прогибов оси по / = 1,2, ..., (п- 1)-Й 


собственным формам, а вызовет лишь дополнительный прогиб оси по / = л-й собственной форме. 
Дополнительный упругий прогиб оси по более высоким формам изгиба при этом будет пренебре- 
жимо малым. 

На рис. 1, а сплошной линией представлена схема исходного искривления оси ротора по 
Г = п-й собственной форме, а пунктирной линией — её искривление по той же форме после при- 


ведения ротора во вращение со скоростью и, <® <. 


Представим исходное искривление оси ротора по л-й собственной форме изгиба как ре- 
зультат её пластической деформации под действием чередующихся по направлению и равных по 
модулю сил Р, приложенных к ротору в л точках на его оси (см. рис. 1, 6). Для описания этого 


искривления нашли реакции Ки = 0,5Р, К, = (-1)”" .0,5Р опор ротора от действия приложенных 


сил. Разбив балку на участки 1, 2, ..., п + 1, записали для каждого участка выражение для изги- 
бающего момента и, использовав известный из [4] метод расчёта деформаций балок, получили 
общие уравнения смещений у, поперечных сечений оси балки, изогнутой под действием прило- 
женной к ней совокупности усилий Р 
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где / = 1...п — номер участка балки, содержащего выпуклую или вогнутую волну л-й собственной 
формы изгиба оси; Е — модуль упругости материала упруго деформируемого вала ротора; 
]— момент инерции сечения вала ротора. 
Экстремальное значение уж„х прогиба оси (стрелы прогиба) на каждой Ёй волне л-й соб- 
ственной формы изгиба оси ротора определили из (1) как 
Р/З 
48Е7? ` 
Из (2) нашли значение сил Р, которыми нагружена балка по схеме рис. 1, 6, которое обес- 
печивает стрелу прогиба уж„х при л-й собственной форме изгиба этой балки: 
‚ 48Е/п? 
в = (-1 [3 У птах г (3) 
Подставив модуль этой силы в (1), получили окончательные выражения для описания ис- 
ходного неупругого искривления оси ротора по л-й собственной форме: 


Уи У 1 фа (-1) (2) 


х=5 +1), 
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Рис. 1. К определению закономерностей упругого изгиба оси ротора по л-й собственной форме 
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Построенная по этим уравнениям линия у, исходного неупругого искривления оси ротора 
по /-й собственной форме изгиба представлена сплошной линией на рис. 1, а. Такое искривление 
оси ротора приводит к возникновению на каждом из его участков длиной //п дисбалансов, рас- 
пределённых с непостоянной интенсивностью: 

ав (х) = ту „ах. (5) 


Используя последнее, нашли значение ДР», главного вектора дисбалансов, обусловленных 
исходным изгибом оси ротора на участке //п. Для этого подставили в (5) первое из уравнений 


(4) и после проведения интегрирования получили 


— Пт __ 1/2п __ 5 т-— 
Ри = | Чы(х) =2т, ] У›аХ = т У птьк, (6) 


Приведение рассматриваемого ротора во вращение со скоростью ®„„,<@<®,„ вызо- 


(1) 
вет дополнительное (к исходному) упругое искривление его оси по л-й собственной форме. Оно 
представлено пунктирной линией на рис. 1, аи вызвано действием лишь инерционных усилий от 
дисбалансов, распределённых по (5) и представленных на рис. 1, в. Найдём этот прогиб, опери- 
руя значениями дисбалансов, а не инерционными усилиями, им пропорциональными. 

Заменив на каждом участке //п распределённый дисбаланс эквивалентным ему главным 


вектором дисбалансов Ди» ‚ определённым по (6), нашли реакции опор: 
Рл=0,5Ри, Бе = (-1)”" 0,55». (7) 


Выделив на левом участке сечение с координатой х! (рис. 1, в), определили равнодейст- 
вующий дисбаланс от распределённой нагрузки при 0 < х<х, — 


РБ», = [тоу „ах я 
о 
что после подстановки первого из уравнений (4) с / = 1 и проведения интегрирования позволило 
получить 
у Топ А- № 
В», [тк т. ' (8) 


Координату плоскости приведения дисбаланса ДР», определили как 


Хх: 
т | 
0 
с =Б- | Х-У@, 
0 


что после подстановки (4), (8) и проведения интегрирования дало 
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2 4 х? 
Хи: =Х, у 2 : (9) 
315 ./ 2 
3—=2Х, 
п 
Изгибающий момент М, в сечении жи с учётом (7)...(9) привели к виду 
т 5 1 п? Е 3. 
М/ = х ХЕ. Х, |. 
и1 Г Ура Е 1 2 [2 1 5 т . (10) 
Дифференциальное уравнение упругой оси ротора на участке 0 < х, < //2п имеет вид 
пет, Р Е 
ах? Е] п птах 16 1 2 [2 1 5 [4 1 
Интегрируя дважды последнее уравнение с начальными условиями: х, =0 у =0, 
ж = 5. > РА = 0, в первом приближении получили уравнение упруго деформированной оси ро- 
п х, 


тора на участке 0 < х, < //2п от действия неуравновешенности, обусловленной исходным неупру- 


гим искривлением этой оси по л-й собственной форме, при приведении ротора во вращение со 
скоростью 9, <<, — 





(1) _ : | . . 
о Г о 


от | А ' 
птах 1 ь 
На участке //2п < х, <//п упругий изгиб оси будет симметричным относительно точки 

х, = //2п , в которой имеет место экстремум упругого прогиба оси (рис. 1, 7: 


2 3 
@°т/ 
а =а, ` Уптах , а, =0,01- +, (12) 
где а, — коэффициент упругого приращения стрелы прогиба оси ротора по л-й собственной фор- 


ме изгиба при его вращении со скоростью и, <® <. 


Последний результат получили подстановкой х =//2п в (11). Анализируя его, отметим, 
что значение а, представляет собой не что иное, как аналитическое выражение квадрата отно- 


шения 6/6 так что 
ЕЛ 
„ =101° Е ЕЙ (13) 


Поэтому для рассматриваемого ротора любая л-я критическая скорость больше первой критиче- 
СКОЙ @нр1 В /7 раз, и поскольку ри -1) < <, ТО а, < 1. 


крп ! 





крп ! 


Как отмечено выше, уравнения (11), (12) позволяют лишь в первом приближении оценить 
дополнительный упругий прогиб оси ротора от исходного её искривления по л-Й собственной 


форме со стрелой у„„„*0 (рис. 1, а). Найденные по этим уравнениям значения у‘) (х,) и 


у® приводят к появлению новых распределённых неравномерно дисбалансов (показаны на 


пу тах 
рис. 1, Г) дополнительно к дисбалансам от исходного неупругого искривления оси ротора по л-й 
собственной форме. Эти дополнительные дисбалансы ведут к новому дополнительному упругому 
прогибу оси ротора по л-й собственной форме и т. д. Учитывая это и используя формулу суммы 
членов бесконечно убывающей геометрической прогрессии, нашли стрелу прогиба оси ротора в 
установившемся его вращении со скоростью в: 
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—(о) — ее к а_- у 
У лутах = Иа, =— пттех , а, < 1. (14) 
К 1- а, 


Заметим, что при а, —>1 стрела дополнительного упругого прогиба оси нежёсткого рото- 


ра по л-й собственной форме у“ >. При отсутствии исходного неупругого искривления оси 


пу тах 


ротора по л-й собственной форме (у„„=0) расчёт по (14) даёт у“ =0 при а, <1. При 


пу тах 
а, =1 тот же случай приводит к неопределённости, отражающей неустойчивое состояние ротора, 


когда бесконечно малое значение у,» приводит к неограниченному росту у“, , характерному 


пу тах ! 
для резонансного режима. 
Поскольку упругий и исходный (неупругий) прогибы оси ротора по каждой л-й собствен- 
ной форме (л = 1, 2, ...) лежат в одной осевой плоскости, суммарная стрела прогиба оси по л-Й 
собственной форме определится как 





(6) (и) у 
Уптах = птах + У лутах = ей 1 
й у у 1 га а, ( 5) 


При а, >1 расчёт по (14) даёт у“, <0 и (15) приводит к У® <у„„. Это согласуется 


пу тах птах 
с теорией колебаний валов, констатирующей уменьшение неупругой деформации валов при их 
вращении со скоростями выше критических. Однако при пуске машины переход и =,„ может 
привести к негативным для всей конструкции машины последствиям и потому недопустим без 
проведения балансировки по л-й собственной форме. 
Для жёстких роторов а, =0, и потому для них у“), =0 и У® =у ть. 


пу тах п мах 
Механико-математическая модель неуравновешенности ротора, балансируемого по 
п > 1-й собственной форме изгиба. На каждом участке ротора, имеющем длину //п, дефор- 


мации оси на скорости ‹& будут приводить к возникновению результирующего локального вектора 
—(«) . 
Ри’ распределённых неравномерно дисбалансов (рис. 1, Г), определяемого с помощью (6) как 


—(& 5 ТЫ) 5т 
Ри = * У пах — ы м У тах 
8 п 





р (16) 


— ©) —(о) 
где У»"ах определено по (15), а а, <1 по (12). Приложен р» в центре каждого рассматриваемо- 
го участка ротора длиной //п (аналогично тому, как приложены векторы Р„ на рис. 1, в). 

р .: —(®) 
Приведём всю совокупность определённых по (16) результирующих дисбалансов Р„ к 


—() —— (6) 
главному вектору Дм и главному моменту М» дисбалансов ротора в целом, имеющих место 











при вращении его со скоростью м, <<: 
п 0+1 
Го] 1 1 Го) [65 1 1 Го] 
р ЕНИБи ме _1-С 9. (17) 
2 2 2 


При принятой на рис. 1 системе отсчёта деформаций оси ротора, когда на первом слева 


а —() 
участке оси протяжённостью //п имеет место у<0, ДР.м будет всегда направлен в сторону 
— - —() 

У пах На этом участке и приложен в центре пролёта / между опорами ротора. При этом Ре =0 и 


—— (6) . 
Мо =0 будет иметь место лишь для нечётных номеров л собственных форм изгиба оси ротора. 


ь — (6) [. уе 
Главный момент М» дисбалансов ротора будет иметь место лишь для чётных значений п при 
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— (=) о ты К К 
Ре" =0. Этот момент лежит в плоскости исходной деформации оси ротора по чётной п-й собст- 


венной форме изгиба и при принятой и оговорённой выше системе отсчёта этих деформаций он 
будет направлен против часовой стрелки. При отсутствии исходного искривления оси ротора по 
п-й собственной форме стрела такого искривления на участке //п оси ротора у„„ =0, и значе- 


(о) о) 
ния Рст, М» обращаются в ноль при любом „<<. 


(п-1) 
Если рассматриваемый ротор при у„„ +0 привести во вращение со скоростью а на ба- 
лансировочном станке, то измерительная система станка зафиксирует в плоскостях коррекции 1, 


2 необходимость установки в них корректирующих дисбалансов Вы, в параллельных при 
нечётном и антипараллельных при чётном п: 
— (в — (в 
58 --5905-С № 0 еоЯ-^, М, ав 
В В В В 
где А, С— расстояния от опор ротора до ближайших плоскостей коррекции; В — расстояние меж- 
ду плоскостями коррекции; /= А+ В+ С. 
Уравнения (17), (18) и соотношения (12), (16), раскрывающие входящие в них величины, 
представляют собой механико-математическую модель неуравновешенности ротора, балансируе- 
мого по (л > 1)-й собственной форме изгиба. Входом этой модели является варьируемое значение 


© <®<0,„. Собственные свойства модели характеризуются значением и углом стрелы 











кр(п-1) 
Усы исходного искривления оси ротора по П-Й собственной форме изгиба, а также параметрами 


т, Е, Л ГА, В, Сконструкции ротора, определяющими по (12) комплексную характеристику а, и 


входящими в (18). Выходы модели — характеристики неуравновешенности ротора о во 


Мы, ри, 5, ), определяемые по (16)...(18) при любом 6 <. 

Идентификация параметров механико-математических моделей неуравновешенности 
роторов и их балансировка. Практическое использование построенной модели неуравнове- 
шенности для решения задач балансировки ротационных агрегатов машин при их проектирова- 
нии, производстве, эксплуатации, ремонте требует знания всех характеристик и параметров соб- 


ственных свойств. К ним относятся отмеченные выше Я — стрелы исходного неупругого про- 


гиба оси ротора по / = 1, 2, ..., п, ... собственным формам изгиба, а также значения 27, Е, 1 / А, В, 
С. Все эти характеристики для полного множества ротационных агрегатов каждого наименования 
являются случайными величинами с различными уровнями рассеивания их значений. 

На проектной стадии создания ротационного агрегата конструктор имеет дело с оценоч- 
ными (практически предельными, наиболее вероятными) значениями характеристик и параметров 
собственных свойств объекта моделирования, что позволяет прогнозировать с помощью модели 
неуравновешенное состояние ротационного агрегата на работающей машине и обеспечивать тем 
самым принятие обоснованных решений по его балансировке. Поэтому все параметры модели 
неуравновешенности проектируемого ротационного агрегата могут быть с той или иной достовер- 
ностью определены по различного рода справочным источникам. При этом в расчётах должны 
приниматься такие их значения, которые обеспечивают прогнозирование практически предельно- 
го возможного неуравновешенного состояния проектируемого ротора. 

Использование модели неуравновешенности для проведения балансировки конкретного 
типоразмера ротора требует знания точных значений всех параметров и характеристик собствен- 
ных свойств модели каждого балансируемого экземпляра ротора. Если этот типоразмер ротора 
необходимо балансировать по / = л-й собственной форме, то он должен проходить балансировку 
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последовательно по / = 1, 2, ..., л-й собственным формам. При этом для каждой из форм должна 
быть построена модель неуравновешенности с точными значениями всех параметров и характе- 
ристик собственных свойств. 

Выражение для а, по (12) можно представить как 


3 
— 72 а т/ — В, 
ВВ ВОО (19) 
получив выражение для комплексной характеристики В„ собственных свойств модели неуравно- 
вешенности ротора, балансируемого на любой л-й собственной форме. По своему физическому 


смыслу В, — величина обратная квадрату &„„ по (13). Эта характеристика определяется пара- 





метрами /т, / Е, Л модели, о которых речь велась выше. При балансировке ротора по любой соб- 
ственной форме точному определению должны подвергаться значения У и В» т. к. именно 
эти характеристики и определяют выходы модели по (16)...(19) как функции <. 

Балансировка по 1-й собственной форме изгиба требует проведения низкочастотной ба- 


лансировки на скорости в, < „„, ‚ что предопределяет специфику балансировки роторов по этой 


форме, теоретически обоснованную и практически апробированную [1]. Там же представлена 
методика идентификации значения В!, входящего в (19). 
С решением последней задачи и расчётом В„ по последнему уравнению, идентификации 


подлежит лишь параметр у„„„ — стрела исходного искривления оси ротора по л-й собственной 


форме изгиба (рис. 1, а). 


Приводя ротор во вращение на балансировочном станке со скоростью &„„, <@ < 


крп ! 


(п-1) 
будем иметь зафиксированными измерительной системой станка значения и углы корректирую- 


—) —(®) р > 
щих дисбалансов Ри:, Рю? в плоскостях коррекции у опор ротора. Для нечётных значений л >1 
эти дисбалансы будут параллельными, для чётных — антипараллельными. 

Имея определёнными эти дисбалансы, исходя из (17), (18), получили общее выражение 


для значения О“) результирующего вектора в} распределённых неравномерно на участке //п 


дисбалансов ротора: 





Си 1+ (-1)’ В 


Б®) — (20 +0) + : 
2 2 1 


ип Кп1 Кп2 


(20) 


Последнее с помощью (16) позволяет идентифицировать значение параметра У. как 


Е ЕЙ 
ты = (2+0) Е: вв | И] (21) 





5 2 2 | п 


Угол у„: вектора у.» на первой волне в системе координат, связанной с ротором, опре- 


деляется как у„, =ф, +п, где ф, — зафиксированный балансировочным станком угол корректи- 


— (о) 
рующего дисбаланса Ди: . 
Последующая балансировка ротора по каждой из / = 2, 3, ..., п форм изгиба сводится к ус- 


тановке во всех л точках экстремума изгиба его оси корректирующих дисбалансов Ро =-Ри, оп- 
ределённых по (6). 


Заключение. Для ротора, имеющего | <, <0,„ и потому требующего последовательной 


(1-1) 
балансировки по / = 1, 2, ..., п-й собственным формам изгиба оси, все параметры его механико- 
математических моделей неуравновешенности при балансировке по каждой собственной форме 
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определяются по представленным выше соотношениям. Эти модели могут успешно использовать- 
ся для решения задач балансировки как при проектировании роторов, так и в их производстве и 
эксплуатации. 
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РУМАМТС $ТМОЕАТТОМ АМО ЕЕЕХТВЕЕ ВОТОВ ВАГАМСТМС ТМ п" ВЕМО $НАРЕ 


О. О. Рои$НКт, О. А. Рои$ВКт 
(Роп За Тесптка! Упмег$Ку) 


5оте теспап/зтз о! аеегттта {те ПехЬ/е гоог ах! еаз8с БисКта т пе едептоде Бепа апа {пегею ге еа 
@зтриеа ипра/апсез аге зресНеа Ёгоидй №5 адупаткс ятиайоп. А пеи тето4 оЁ еп Иу/тпд {пе тодЕ! рага- 
тег апа Бейапстд оЁзисй гоюг$ п {пе едептосе ВБепа [5 емо/ед. 

Кеуигогаб5: ПехЬ гоогз, дупати/с$, Баапста, тодЕ/тд. 
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